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Contributions of the direct quark mechanism are studied in the nonmesonic weak decays of 
hght hypernuclei. The AA^ NN transition is described by the one pion exchange mechanism 
and the direct quark mechanism, induced by the four-quark vertices in the effective weak Hamil- 
tonian. By employing a realistic wave function of the A inside the hypernuclei, nonmesonic decay 
rates of y^^H, jjie, and ^He are calculated. The results show that the direct quark mechanism is 
significantly large and gives a large AI = 3/2 contribution in the J = channel. The relative 
phase between the one-pion exchange and the direct-quark contributions is determined so that the 
effective weak Hamiltonian for quarks give both of them consistently. We find that the sum of these 
two contributions reproduce the current available experimental data fairly well. 
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1. Introduction 



An 5 = — 1 hypernucleus decays from its ground state to non-strange hadrons, emitting roughly 
170 MeV (~ A/a — Mat) of extra energy[l|. As is the decay of free A, it can emit a pion + 30 MeV. It 
is, however, known that the pionic decay is strongly suppressed because the remaining energy is too 
small for the produced nucleon to go above the Fermi surface of the nucleus. Consequently, the main 
decay mode of (heavy) hypernuclei is the decay without a pion emission, called nonmesonic decay. 
The nonmesonic decay releases the full 170 MeV of energy and therefore is not Pauli blocked. The 
final state is dominated by two nucleons ejected in the opposite directions so that the momentum 
is conserved, although three-body processes may not be negligible [Q]. A natural way of describing 
the nonmesonic decay is to assume that the emitted (virtual) pion from the A decay is absorbed 
by a nucleon in the nucleus. This is called the one-pion exchange (OPE) weak transition and has 
been studied by many authors [Q,^^,^^,^]. 

On the other hand, the final A^A^ state has a large relative momentum so that the short-range 
interactions are also important ||, 10 1. Two of the present authors have proposed that the quark 
structure of the baryons provides a new mechanism for the nonmesonic AA^ NN, called direct 
quark (DQ) mechanism |ll , 1 2 , 13| Jl^ . We calculated transition amplitudes for the S wave initial 
states in both OPE and DQ, and found that the contribution of the DQ is in general as large 
as that of the OPE. DQ is even dominant in the L = to L = transitions. It gives a large 
contribution in the J=0 (L=0,S=0) transitions and thus contributes to the nn (1=1) final state. 
We found that the nn/pn ratio of the final state becomes larger than that predicted in OPE. It 
was also found that the J = transition amplitudes have a significant AI = 3/2 contribution. 

In this paper, we apply the DQ transition potential to the decay of three light hypernuclei, \H., 
^He, and ^He in order to study how the features of DQ are reflected in the decay observables of 
simple hypernuclei. We employ realistic wave functions for the hypernuclei, and evaluate the matrix 
elements of the two-body weak transition potentials in OPE and DQ, and also of the potential given 
by the superposition of OPE and DQ. The relative phase of OPE and DQ are determined so that 
the effective weak Hamiltonian for quarks give both consistently. The results are compared with 
current available data. We find that the OPE+DQ superposed transition potential reproduce the 
currently available experimental data fairly well. 

In sec. 2, we present the two-body weak transition potentials for AA^ — > A^A^ in DQ and OPE. 
In sec. 3, the nuclear wave functions are given. In sec. 4, the results of our calculation are presented 
and are compared with experiment. Sec. 5 is devoted to conclusions. 
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Figure 1: The quark diagrams for the AN NN transition. • denotes a weak 4-quark vertex. 
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Figure 2: The quark diagrams for the AA^ NN transition. • denotes a weak 4-quark vertex. 



2. Weak transition potential for AA^ —>■ NN 



Considering that a baryon has three constituent quarks, AA^ NN process can be described 
by the diagram, such as those shown in Fig. || and Fig. ^. In the diagrams shown in Fig. |l|, 
the strangeness changing weak interaction and an emission of a constituent quark anti-quark pair 
take place in the A hyperon, and the pair is absorbed by a constituent quark in the nucleon. The 
diagrams in Fig. ^ show the weak interaction of two constituent quarks in a totally anti-symmetric 
six constituent quark state. The first picture might be well represented by the diagram in Fig. ^ 
where the baryon is a Dirac particle and couples to, for example, a pion by a phenomenological 
Yukawa type vertex. This is the one called the meson exchange mechanism and has been studied 
welip,^j5|,^j^,p. Though this picture is very natural, one sees that this picture cannot be valid in 
the region where the two baryons overlap with each other. In such a region, the diagrams in Fig. 
^, namely the direct quark process, might cause the AA^ NN transition. 

Recently, Inoue, Takeuchi and Oka pointed out the importance of the direct quark (DQ) pro- 
cesses and derived the DQ transition potential for the AA^ NN systems [|ll, 12 There we 
employed the nonrelativistic constituent quark model for the baryons and the two baryon states are 
constructed according to the quark cluster model, which takes the full quark antisymmetrization 
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Figure 3: One pion exchange diagram. 

into account, 

\BiB2ik, L, S, J)) = A\cf>icf>2Xik)iL, S, J)) (1) 

where (f)i and (p2 denote the internal wave functions of the baryon Bi and B2 and x{k) describes 
the radial part of the relative motion of the baryons with k relative momentum. A is the quark 
antisymmetrization operator. 

The weak interaction is represented by the effective Hamiltonian, 

^eir=-^ E KrOr, (2) 

where 





K2 


K3 


K5 


Kq 


-0.284 


0.009 


0.026 


0.004 


-0.021 



and 

01 = {daSa)v-A{ui3Ui3)v-A - iUaSa)v-A{dpUp)v-A (3) 

02 = {daSa)v-Aiui3Ui3)v-A + {UaSa)v-Aidl3Ui3)v-A 

+2{daSa)v-A{d/3di3)v-A + 2{daSa)v-Ais/3S/3)v-A (4) 

03 = 2{daSa)v-A{ui3Ui3)v-A + 2{UaSa)v-A{dl3Ui3)v-A 

-{daSa)v-A{dpdp)v-A - {daSa)v-A{spSp)v-A (5) 
O5 = {daSa)v-A{ui3Ui3 + J/jd/j + (6) 
Oq = {daSf})v-A{upUa + dpda + SpSa)v+A ■ (7) 
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This is derived from the W exchange diagrams in the standard theory, by taking perturbative 



corrections due to the strong interaction] 17]. We evaluated the two-baryon matrix elements of the 
effective weak Hamiltonian H^j^^^, using the six-quark wave functions given in eq.(|l|), 

Vik,k') L,,s,,j ={NN{k',Lf,Sf,J)\H^ff'\ANik,U,Si,J)) . (8) 

Lj,Sj„J 

Here, k, L, S, J stand for the relative momentum of two baryons, the orbital angular momentum, 
the spin, and the total angular momentum of the initial and the final states, respectively. The 
obtained two-body transition amplitudes are regarded as the transition potential in the momentum 
representation for each channel specified by L, S and J. Fig. ^ shows various diagrams contributing 
to the transition potential, where the dots represent the weak four quark vertices. The explicit forms 



of the transition potential is given in ref.[14|. 



In the present study, we calculate the decays of the s-shell hypernuclei from its ground states. 
If we employ the simplest shell model wave functions for the hypernuclei, we only need the L, = 
transition potentials. The relevant transition channels are given in Table 1. The direct quark 
induced transition potential depends on two quark model parameters that are the constituent 
quark mass, m, and the Gaussian parameter, b. We use m = 313 = M]\f/3 MeV and b = 0.5 fm in 
the present calculation. 

As for the meson exchange weak transition, we here consider only the one-pion exchange (OPE) 
process, and will discuss other contributions later. The pion exchange transition potential is ob- 
tained by evaluating the diagram like Fig. |3|. We here employ the following strong and weak pion 
vertices, 

Han-k = i9wi^n{^ + A75)7r°VA - iV2gy,i)p{l + A75)7r+VA (10) 

where the strong coupling constant Qs is chosen as the standard value of the ttNN coupling: 
Qs = —13.26. The weak coupling constants Qw and A are determined so as to reproduce the 
A — > Nt: decay amplitudes: 

= -2.3 X 10^^ (11) 
A = -6.9 . (12) 

The relative sign of gs and gyj is important when we consider the interference of OPE and 
DQ contributions. Namely, gyj must be chosen consistently with our weak Hamiltonian for quarks 



and the baryon wave functions. In order to check the consistency, we use the soft pion theorem. 
Consider the parity- violating part of the A — > mr^ decay matrix element in the soft pion limit, 

lim{n7r%q)\H,ff{PV)\A) = -f (n| [QlH.jfiPV)] |A) . (13) 

The soft-pion limit is taken by assuming 

(^O(g)l A«(:e) |0) = -if^q^expiiq ■ x) . (14) 

Note that the Goldberger Treiman relation 

Jtt = 7 7 (15) 

is satisfied, where Qs is the strong ttNN coupling constant defined in eq.(^). Because Hf,ff{AS = 1) 
consists only of the left-handed currents and the flavor singlet right-handed currents, we have 

[Ql H,ff{PV)] = [Ql - Ql H,ff{PV)] (16) 

= -[Q% + QlHeff{PC)] = -[l3,Heff{PC)] (17) 



Thus we obtain 



lim (n^O(g)l H^ffiPV) |A) = H,ff(PC) |A) (18) 



Comparing this with the effective hamiltonian eq.(lO), the weak ttNN coupling constant is given 
by 

g^ = ^^{n\H,ff{PC)\A) <0 (19) 

after evaluating the matrix element in the right hand side, which is positive with our definition of 
the Hamiltonian H^^^j^^ and the quark model wave function. 

In Table we summarize the radial part of the OPE transition potential for the nine possible 
channels which appear in the calculation for the s-shell hypernuclei. The function in the table are 
defined by 

r/ \ 9w9s l^n ~ 6 ^ /onN 
fv) = — ; ^""T-TT (20) 

Vir) = 1 + ^ (21) 

.sill 
T r = - + ^ + 22 

6 rriTrr mtr'^ 
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Table 1: One pion exchange induced transition potential 

spin-orbital /jv Potential 

~ap pk-^pn ^So ^So 1 --^A^/(r)~" 

hp ^'Po 1 i^V{r)f{r) 

dp -^'D, _6A^r(r)/(r) 

ep ^iPi iy|y(r)/(r) 

On 



3. 

nA^nn ^Sq^^Sq 1 -A^/(r) 



2M 

bn ^^Po 1 i^O/W 

/n 'S.^'P, 1 _iy|y(r)/(r) 



when the form factor is not taken account. Here M is the average mass of the baryons and 



TTiTr = Y m'^ — Qq is an effective pion mass introduced in order to take care of the finite energy 
transfer. The radial functions are modified when the form factor with a cutoff A,r is introduced, so 
that 

/W-/(r)- f4^)'/(A.r) (23) 
f{r)V{r) ^ f{r)V{r) - (^) f{A^r)V{A^r) (24) 
f{r)T{r)^ f{r)T{r)-{^\ f{A^r)T{A^r) . (25) 

In the present calculation, we choose A| = 20m'^, according to ref.[^. 

We regard DQ and OPE as independent of each other. Fig. ^ shows the various diagrams which 
contribute to the DQ transition potential. Among them, the diagrams (b) - (g) contain the quark 
exchange between the two baryons, and therefore the range of the DQ potential is determined by 
the size of the quark wave function of the baryon. Thus, their contributions are independent from 
the one-pion exchange, as the strong repulsion due to the quark-exchange is independent from the 
one-pion exchange potential in the nuclear force. 

As the diagram (a), on the other hand, gives a long-range contribution, it may cause a double 
counting problem. This is the process in which A decays into n with another nucleon as a spectator. 



B3 B4 B3 B4 




Bi B2 Bi B: 

Figure 4: The direct quark transitions. • denotes a weak 4-quark vertex. 

Its matrix element is nonzero only through high momentum components of the two baryon wave 
functions and therefore comes mainly from the short-range correlation of the baryons. Because 
the one-pion exchange seems not responsible for the short-range correlation, we also regards the 
diagram (a) as independent. It should also be noted that the contribution of the diagram (a) is in 
general small as is expected. 

3. Initial and final wave functions 

For the s-shell hypernuclei, we employ a shell model wave function for the nucleon core, denoted 
by cluster, x. It has been pointed out that it is important to solve the A single particle state carefully 
using the realistic A- nucleon interaction p^]. We denote the relative motion of the A hyperon against 
the nucleon cluster x by a wave function ^{rjC) with ta, relative coordinate. The Hamiltonian H 
of the K-x system is given by the sum of the kinetic energy T and a potential energy V\x- The 
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potential energy V\x is obtained by folding the potential energy between A and the nucleons. In 
the present study, we construct it by folding the YNG potential with (Os)^ where x stands for the 
number of nucleon in the cluster x. The YNG potential is a three-range Gaussian potential which 



reproduces the G- matrix for the Nijmegen AA^ potential model D|24]. The Fermi momentum kf is 
fixed 0.9 fm^^ for light hypernuclei. 

The Schrodinger equation is solved variationally using the local Gaussian basis functions. We 
expand the wave function for the state which has orbital angular momentum / and the spin S as 

^{rA) = ^fi,id)Ml;d),S;j) . (26) 
d 

where j = I + S. The spin of hypernuclei, S, is the sum of total angular momentum of the cluster 
X, Jx, and the spin of A, S\. In eg. ([2^), <!>(/; d) is the basis function given by 



^l;d)=UrA;d)Yi{rA) (27) 
MrA;d) = 47r(V^bA,.)"^/^exp[ ~^l+^'V /(^) (28) 



26ax '^''6^ 



Ax 



The parameter d can be considered as a generator coordinate describing the distance between two 
clusters. The Ji is the l-ih order modified spherical Bessel function. The choice of b\x as in 
eq(p9|) enables us to use the Talmi-Moshinsky transformation coefficients in calculating the two- 
body matrix elements later. The variational amplitudes fij satisfy the generator coordinate method 
equation 

[Hij {di , ds) - ENij {di , d2)] fij {d2)=0 (30) 

d2 

where the energy and the normalization kernels are given by 

Hij{di,d2) = m;di),S;j\HMl;d2),S;j) (31) 
Nij{di,d2) = m;di),S;j\lMl;d2),S;j) . (32) 



For the hypernucleus we take the Gaussian parameter for Os as 1.358 fm[25,2f:]. We 



present the folding potential for even / state in Fig |5[ With the following seven values for d, 

d = 0.0, 1.0, 2.5, 4.0, 5.5, 7.0 and 8.0 fm. (33) 
we obtain the ground state energy E = —3.08 MeV and corresponding variational amplitudes. 



V [MeV] 



-10 
-20 




Figure 5: Folded potential for ^He for even / state. 

Fig 1^ shows the obtained wave function ^{r). This wave function is considerably different 
from the one which is obtained with the one range Gaussian (ORG) potential. As shown in Fig |5|, 
the potential V\x is repulsive at the short distance, reflecting the presence of a repulsive soft core 
in the YNG interaction. Therefore the obtained wave function is pushed out. This feature gives 
rise to a sizable difference in the r.m.s. distance between a and A, 3.06 fm for YNG while 2.69 
fm for ORG 1 25, 26 1 . Motoba et al. applied this wave function to the pionic decay of the s-shell 
hypernuclei and found that the YNG wave function gives a better agreement to experiment than 



the ORG one 1 27 



The wave function ^[r\) for ^He (or ^H) is calculated in the same way. In the 4-body case, 



the Gaussian parameter fe^r is taken as 1.65 fm[25,26|. 



The AA^ wave function, ^'aat, is given by the product of the above wave function, '^{r\), and 
nucleon wave function, ^Ar(rAr), which is assumed as a Gaussian, 

^kN{rK,rN) = ^{rx) x ^N{rN) ■ (34) 

We perform the Talmi-Moshinsky transform for the two-baryon AA^ wave function ^ \n in order 
to separate the relative motion and the center of mass motion. Defining r = r\ — and the 
two-baryon center-of mass coordinate R, the product wave function can be expanded as 

'$AN{r,R) =J2WnNUno{r,br)Yoo{r)um{R,bR)Yoo{R) , (35) 



nN 
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Figure 6: Wave function for A-a system, 
where Uno is the harmonic oscillator eigenfunctions for {n,l = 0), and 



^ = J^^^^±^6^ (36) 



Ma + xMn , 

The b parameters, A=4 are 6,- = 2.24 fm, bn = 1.61 fm, and for A=5 are br = 1.84 fm, b^ = 1.21 
fm, respectively. 

We consider the short-range correlation between A and N. In evaluating the OPE transition 
amplitude, we use the short range correlation obtained for the Nijmegen A-A^ potential model 



D|25j2q|. The Nijmegen model D, however, has a hard repulsive core at short distance, which 
seems inconsistent with the constituent quark picture. Therefore, in evaluating the matrix element 
of the DQ transition potential, we use the following form of the short range correlation. 



1 / , 

Uno{r, br) j^Uno{r, br) \l - Ci exp[--2] ) . (38) 



Here, Ci and rj are parameters, and is a normalization factor. 

As wave function of the outgoing two nucleons in the final state, we employ the scattering 
state obtained in the Nijmegen model D[p^,p^]. On the other hand, we use the plane wave with 
the following form of the short range correlation in evaluating the DQ contribution because of the 
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same reason stated above. 

.2 

r 



Mkr) ^ jiikr) ( 1 - C/exp[-^]j (39) 



We evaluate the DQ matrix elements for several choices of the parameters for the short-range 
correlation, while we assume Cf = Ci and rj = ri for simplicity. 

4. Results and Discussion 

In this section, we present our results and compare them with experiment. The experimental 
data are taken from ref ^,|0|,|l|,|2|,||,|4|,||,||] . 



The nonmesonic decay rate of the hypernucleus is given by 

r = ^<5(^.(:7.)|(*/|i/""^"'=|*i)|2 (40) 
/ 

where 6{E.C.) stands for the delta function for the energy conservation. For the final state, we 
observe only the two nucleons outgoing in the decay. Denoting the momenta of the nucleons as ki 
and /c2, we find 

J {2-kY (27r)^ tl^ 



and 



5{E.C) = 6 (^Mres + 2Mm + ^ + ^ - Mres - Ma - Ea - Mn - En^ . (42) 
Here V^'^"'^ is a two-body transition potential, and ^'j can be expressed as 

= E ^res X ^AN{ch) ^ (cfp for initial state for channel ch) (43) 

ch 

for which the spin is implicit. It is easy to obtain the (cfp) for the initial state in the present 
(Os)^' -|- A configuration, resulting a factor for the number of pairs and their spin average factor for 
each channel, labeled by a-f in Table 1, of the AA^ states. It should be noted that these channels 
do not interfere with others. Thus we obtain F = J2ch^ch, and the decay rate for each channel is 
given by 

Tc/i = of AA^ pairs] x [spin average factor] x (44) 
"^^^ '^^\27T6iE.C.)\{ki,k2\V,T''''\AN)\^ . (45) 



(2^)3 (27r)3' 
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4.1. ^He 

In ^He, there are two Kp bonds and two An bonds. The Kp or An system takes both spin S=0 
and spin S=l state. Thus the spin average factor becomes 1/4 (3/4) for S = Q {S = \) channel. 
The total nonmesonic decay rate of ^He is given by 

r„^(iHe) = + + + + + + + + (46) 

where T^ ~ T^ are the partial decay rates and calculated by 

^ a,b 

c,d,e,/ 
a,b 

First we study the OPE mechanism. Table |2| shows the partial decay rates obtained for the one 
pion exchange potential. All the decay rates are written in the unit of Fa, the free A decay rate. 
We list three sets of results. The values listed under "OFF" are the results when we omit the form 
factor and the short range correlation into account. The values listed under "FF" are the result 
when we take only the form factor into account. The values listed under "FF and SRC" are the 
final result with both the form factor and the short range correlation taken into account. Because 
A/ = 1/2 is assumed for the weak AA^vr vertex, these partial decay rates satisfy the isospin relation 

^n/^p — bfi/bp — /n/ fp — 2. 

In the "OFF" and "FF" cases, the channel dp is dominant. This comes from the tensor part of 
the transition potential, which is enhanced due to the large relative momentum in the final state. 
One can see that the form factor reduces most of the partial decay rates significantly. In the "FF 
and SRC" case, the channel Cp becomes large and is dominant. This comes also from the tensor 
part of the transition potential together with the tensor part of the final state interaction. In the 
OPE mechanism, the channel Cp has also a large rate, while the rates in the channels a and b are 
very small. 

Table |2| shows calculated nonmesonic decay rates of ^He in OPE. We list the proton induced 
decay rate (Fp), the neutron induced decay rate (r„), the total decay rate {Tnm = Tp + Tn), and 
the n-p ratio {Rnp = r„/Fp). The experimental data are listed in Table ^. One can see that the 
calculated Tp is in good agreement with experiment, while the calculated F„ is much smaller than 

13 



1 

^4 



d^k 
(2^ 

d^k 
(2^ 

d^k 
(2^ 

d^k 



47 (27r)3 



^(27,)6{E.C.)\{NN\Vl,\AN)\' 
^^{27r)6{E.C.)\{NN\V^^^,^^JAN) 
^(27,)6{E.C.)\{NN\V:^,\AN)' 
^{27:)6{E.C.)\{NN\Vp\AN) 



(47) 
(48) 
(49) 
(50) 



Table 2: Calculated nonmesonic decay rates of ^He in the one pion exchange mechanism (in the 
unit of Fa, the decay rate of the free A). 



Ch 




1^ 1^ 


FF and SRC 


Up 


r\ f\r\f\ A 

0.0004 


0.0101 


0.0002 


bp 


O.OlZD 


O.OObO 


0.0031 


Cp 


0.0013 


0.0303 


0.1022 


dp 


o.oyio 


0.1789 


0.0415 


€p 


U.iloz 


U. 1)544 


0.0346 


J" 

Jp 


O.Ozol 


U.Uizi 


0.0093 


On 


U.UUUy 


U.UZUZ 


U.UUUo 


bn 


0.0251 


0.0121 


0.0063 


fn 


0.0502 


0.0242 


0.0186 




0.544 


0.291 


0.191 




0.076 


0.056 


0.025 


r 

nm 


0.621 


0.348 


0.216 




0.140 


0.193 


0.130 



experiment. The Fp is dominated by large contribution of the channels c and d, which vanish in 
the neutron induced decay. Thus the calculated n-p ratio is much smaller than the experimental 
one. 

We turn to the DQ mechanism. Table ^ shows the partial decay rates when we employ the direct 
quark induced potential. We list four sets of result which are calculated with different short range 
correlations, while we assume ri = rj and Ci = Cf = C for simplicity. The parameter tq is fixed as 
0.5 fm, which is equal to the Gaussian parameter b in the quark model. The partial decay rate for 
channel dp is zero, because the direct quark potential has no AL = 2 part due to the nonrelativistic 
truncation at the first order in p/m. The total decay rates, Fp F„ and Tnm decrease as the short 
range correlation strength Ci or Cf increases. This is natural because the short range correlation 
reduces the inner part of the wave functions where the direct quark mechanism is important. 

Compared to the results in Table 2, we find that the overall magnitudes of the decay rates are 
comparable to OPE. The components, however, are different. For the proton induced decays, OPE 
is dominated by the / = final states, Cp, dp, and Cp, while the DQ gives a large decay rates to 
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Table 3: Calculated nonmesonic decay rates of ^He in the direct quark mechanism (in the unit of 
T\, the decay rate of the free A). 



Ch 


C = 


C = 0.3 


C = 0.5 


6 = 0.7 


Op 


0.0130 


0.0150 


0.0167 


0.0185 


bp 


0.0127 


r» m on 

0.0120 


n n 1 1 o 

0.0113 


n ni n cr 

0.0105 


Cp 


0.0968 


0.0690 


0.0548 


0.0440 


dp 














Bp 


O.OOoD 


0.0061 


0.0064 


n r\r\r"~7 

0.0067 


r 

Jp 


0.0345 


n no o 

0.0353 


n no c o 

0.0353 


n no c! o 

0.0352 


an 


U.U i Li 


U.UuiD 


U.U4U / 


n riQoo 
U.OoZZ 


bn 


0.0059 


0.0066 


0.0069 


0.0073 


fn 


0.0622 


0.0642 


0.0648 


0.0650 




0.163 


0.137 


0.125 


0.115 


r„ 


0.141 


0.122 


0.112 


0.104 


r 

nm 


0.304 


0.260 


0.237 


0.219 




0.865 


0.889 


0.903 


0.910 



the 1 = 1 final states, Op, 6p, and fp. The neutron induced decays go only to the 1 = 1 states, a„, 
bn, and /„, and thus the DQ gives much larger contribution than OPE. It should also be stressed 
that the J = channels play major roles in DQ, while the OPE is dominated by J = 1. It is 
important to note that these general features can be tested in the decays of the s-shell hypernuclei. 
For instance, we will see that DQ must be dominant in the decay of because the / = final 
states are prohibited there. 

Next we investigate the AI property of the direct quark mechanism. In the following, we employ 
Ci = Cf = C = 0.5. Table ^ shows the effects of the AI = 3/2 part of the potential. The values 
listed under "DQ full" are the results when we employ the full direct quark induced potential. 
While the values listed under "A/ = 1/2 only" are the results when we omit the AI = 3/2 part. 
One sees that the channels a and 6, which are J = transitions, get significant contribution from 
the AI = 3/2 part. This indicates that AI = 1/2 rule is broken significantly in the non-mesonic 
decays of light hypernuclei. 

The large AI = 3/2 transition is very interesting. It is, in fact, expected naturally in the 
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Table 4: Calculated nonmesonic decay rates of ^He in the direct quark mechanism (Cj = Cf = 
C = 0.5). The partial rates, Cp, dp and Cp do not have AI = 3/2 contribution because the final 
state has 1 = 0. 



Ch 


AI =1/2 only 


DQ full 


Up 


0.0118 


0.0167 


bp 


0.0001 


0.0113 


fp 


0.0334 


0.0353 


On 


0.0237 


0.0407 


bn 


0.0003 


0.0069 


fn 


0.0668 


0.0648 



standard theory of the weak interaction. But one may wonder whether the present quark model 
with Hamiltonian H^^^j^^ is capable to reproduce the AI = 1/2 dominance of the free hyperon 



decay. In order to check it, we study the free A decay with this effective Hamiltonian[37|. We find 
the dominance of quark diagrams with internal weak vertex in which AI =3/2 amplitudes vanish, 
according to the Pati-Woo theorem which is resulted from the color symmetry of the ground state 
baryonsp8|. Therefore, the dominant decay of the hyperons go through the AI = 1/2 part of the 
Hamiltonian. It is extremely interesting to confirm the strong AI = 3/2 weak transition. 

The calculated Fp, F^, Tnm, and Rnp in DQ are shown in Table ^ The proton induced decay 
rate Fp in DQ is smaller than that of OPE. On the other hand, the neutron induced decay rate F„ 
in DQ is much larger than that of the OPE. This is due to large contribution of channel a„ and /„. 
The total nonmesonic decay rate in DQ is roughly equal, while the n-p ratio in DQ is much larger 
than that of OPE and is about 0.9, which is closer to the experimental data. 

As was argued in the section 2, the final results are given by the superposition of the OPE and 
DQ processes. Table | summarizes the results given by the sum of the two transition potentials. 
We find that the neutron induced decay rate, F„, is significantly enhanced from OPE and becomes 
consistent with experimental data. On the other hand, the combined result overestimates the 
proton induced decay rate, Fp. Thus the total decay rate, Tnm-, is slightly overestimated. The 
n-p ratio, Rnp, is predicted in between the values for OPE and DQ, while the experimental data 
suggest larger values. 
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Table 5: Calculated nonmesonic decay rates of ^He (in the unit of Fa). 



Cn 


OPE only 


DQ only 


OPE + DQ 


EXP 


2£] 


EXPpOl 


ap 


0.0002 


0.0167 


0.0188 








7 

bp 


0.0031 


0.0113 


0.0026 








Cp 


0.1022 


0.0548 


0.2612 








dp 


0.0415 





0.0415 








6p 


0.0346 


0.0064 


0.0207 








fp 


0.0093 


0.0353 


0.0763 










0.0003 


0.0407 


0.0356 








I, 

On 


U.UUOo 


U.UUDy 


U.UZD4 








fn 


0.0185 


0.0648 


0.1437 










0.191 


0.125 


0.421 


0.21 ± 0.07 


0.17 ± 0.04 




0.025 


0.112 


0.206 


0.20 ± 0.11 


0.33 ± 0.04 


r 

J- nm 


0.216 


0.237 


0.627 


0.41 ± 0.14 


0.50 ± 0.07 




0.132 


0.903 


0.489 


0.93 ± 0.55 


1.97 ± 0.67 



4.2. ^He 



In the ^He, there are two Ap bonds and one An bond. The An pair is in the spin S=0 state so 
that the spin of ^He is equal to zero. Thus there is no contribution from the channel /„. The the 
partial decay rates of are given by 



3 r d^k r d?K , I, ,„ ,.,2 

(52) 



j^jj^{2^nE.C.)\{NN\V:^,\kN)\ . (53) 

Table ^ shows our results and the experimental data. The results of OPE are qualitatively the 
same to \iie, but are reduced to 80 % or less. On the other hand, the results of DQ have qualitative 
differences. The rates in the channels Cp and a„ are very small, which are large in the ^He case. 
This indicates that the DQ contribution is sensitive to the wave function of the initial AA^ system. 
Again the full calculation of OPE+DQ provides a good agreement to the experimental data except 
that Tp is again slightly overestimated. In our result, r„ is dominated by the channel b. 
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Table 6: Calculated nonmesonic decay rates of ^He (in the unit of r^). 



Cn 


OPE only 


DQ only 


OPE+DQ 


EXP 


35 


1 

1 


Op 


0.0001 


0.0183 


0.0214 








7 

bp 


0.0021 


0.0087 


0.0023 








Cp 


0.0818 


0.0004 


0.0884 








dp 


0.0321 





0.0321 








Gp 


0.0224 


0.0048 


0.0140 








fp 


O.OUoo 


0.0275 


0.0567 








an 


0.0005 


0.0013 


0.0004 








u 

On 


U.UUoo 


U.UiUo 


U.UooU 








fn 



















0.145 


0.060 


0.214 


0.15 ± 0.02 


r„ 


0.009 


0.012 


0.038 


0.04 ± 0.02 


r 

^ nm 


0.154 


0.072 


0.253 


0.19 ± 0.04 




0.061 


0.202 


0.178 


0.27 ± 0.14 



4.3. iH 



In ^H, there are one Ap bond and two An bonds. The Ap pair is in the spin S=0 state so that 
the spin of ^He is equal to zero. Thus there are no contributions from the channels Cp though fp. 
The the partial decay rates are given by 



d^k r d^K 



2 



"^-^ = 7(^7 (2^('")^(^-^-^ pA^KJAiV) (54) 
= 27 / 777^ / T^T^CimE.C.) \{NN\V:AAN) (55) 



47 (27r)3 7 (27r)3 

^/ =^\jwfl Y^m\m\" ■ (56) 

Table ^ shows our results. In our model, each partial decay rate is the same as that of ^He 
except for the number of bonds and the spin average factor. The effect of the interference of OPE 
and DQ is large in the channel /„, which does not appear in the ^He decay. At present, the 
experimental data are very limited. Only the total nonmesonic decay rate is barely known. Our 
prediction agrees with it. The result also suggests that the neutron induced decay is stronger than 
the proton. Thus we expect a large n-p ratio for ^H. We anticipate new good experimental data 
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Table 7: Calculated nonmesonic decay rates of (in the unit of r^). 



Ch 



OPE only 



DQ only 



OPE+DQ 



EXP|6| 



0.0003 
0.0041 



0.0366 
0.0174 



0.0429 
0.0046 



On 
bn 
fn 



0.0003 
0.0041 
0.0130 



0.0006 
0.0054 
0.0504 



0.0002 
0.0190 
0.0107 



nn 



R 



np 



0.004 
0.017 
0.022 
3.952 



0.054 
0.057 
0.110 
1.048 



0.047 
0.126 
0.174 
2.660 



0.15 ± 0.13 



for the nonmesonic decay of ^H. 

4.4. Other contributions 

So far we considered only the one-pion exchange and the direct quark transitions. Several 
other possibilities are considered in the following. First we discuss the effect of the heavy meson 
exchanges. Ramos et al. calculated the AA^ — > A^A^ transition in a full one-boson-exchange 
mechanism [l^] . They include the other pseudosclar mesons, r/ and K, as well as the vector mesons, 
p, 00 and K* . In constructing the transition potential induced by the exchange of the heavy mesons, 
the Nijmegen or the Jiilich strong vertices, Hnnt], Hnnk, Hjyp, H^j^f^^ and Hnnk*-: are employed. 
On the other hand, the weak vertices, H\j\fri, H\j\fK, H\Np, H\j\fi^ and H\nk*-, which cannot be 
determined from the hyperon-decay experiment, are determined by the SU{6) symmetry and the 
soft meson theorem for the PV vertices and pole model for the PC vertices ||20|]. The nonmesonic 
weak decay of j^C is studied in the shell model framework. The result shows that the combined 
TT + p exchange predicts a very similar total decay rate to the OPE one. It is found that the total 
rate is reduced by more than 40% when K exchange is added. It is also found that K* exchange 
compensates the K exchange and that the total rate for the combined n + p + K + K* is only 10 
% smaller than the OPE only. The contribution of r] and lo are very small and tend to cancel each 
other. Therefore the total rate for the combined n + p + K + K* + r] + uj is similar to the OPE only. 
They found that the addition of K exchange reduces the n-p ratio considerably while the addition 
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of other mesons does not change the result much. The final n-p ratio is smaller than that of OPE, 
which is much smaller than the experimental data. 

Shmatikov studied the 27r exchange contribution. It is indicated that the diagrams with TiN 
and A^A^ intermediate states cancel each other and the net effect contributes only to the J = 
amplitudes 1 21]. 



Recently, Itonaga et al. [26| have studied the 2tt/p and 27r/(T exchange weak potentials in which 
27r are coupled to p and a, respectively, in the exchange process. The weak vertices H\n-k and 
Hy.n-k are empirically known from the pionic decay of A and S hyperons and the Lee-Sugawara 
relations. The strong coupling constants are determined so that the corresponding strong potential 
versions of the 2t: / p and 2t: /a exchanges can simulate the OBE potentials due to p and a exchange, 
respectively. There is no phase ambiguity between OPE and 27r/p and 2ti/(t exchange weak poten- 
tials. The 27r/p exchange potential has a tensor part whose sign is opposite to that of the OPE one 
and thus it tends to weaken the tensor interaction of the latter potential. The 27r/(7 exchange is of 
central type typical from its character. These potentials are applied to the non-mesonic decays of 
typical hypernuclei raging from s-shell to medium-heavy systems. The addition of 27r/p exchange 
decreases the non-mesonic decay rates by 5-10 % with respect to the OPE estimates. The inclusion 
of both 2t[ / p and 2T:/a exchanges leads to increase the decay rates of about 5-15 % with respect 
to the OPE estimate. Calculated non-mesonic decay rates, together with minor contributions from 
the pionic decay channels, are generally consistent with the lifetimes measured for p- and sd-shell 
hypernuclei. As for the n/p ratios, the addition of the 27: / p and 2ti/(t exchanges does not cause 
much improvement and the theoretical values are still far from the experimental data. 

In all, the meson exchange contributions other than OPE seem to be less important than DQ. 
Therefore describing the AA^ — > A^A^ transition by the sum of OPE and DQ, seems to be reasonable. 



5. Conclusion 

It is shown that the one pion exchange mechanism is significant to the nonmesonic decay of light 
s-shell hypernuclei but is not able to reproduce some of the experimental data. Our result shows 
that the direct quark processes in AA^ — > A^A^ provides as large contribution as the OPE. We also 
show that the DQ contribution is qualitatively different from the one pion exchange mechanism. 
The DQ contribution is dominant in some channels, such as the proton induced J = channels, 
ap and bp. It is also shown that the direct quark mechanism causes a large A/ = 3/2 transition 
in several channels. Our results are qualitatively consistent with those of Maltman and Shmatikov 
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[10], although the calculated amplitudes have quantitative differences. 

In this paper, we have determined the relative phase of OPE and DQ by using the soft pion 
relation. Thus we superpose these two mechanisms without ambiguity. We reproduce the present 
data for ^He, ^He and fairly well by superposing OPE and DQ. Some predictions are made for 
the nonmesonic decays of ^H, to which experimental data are very limited. It should be stressed 
that the ratios of the partial rates for various channels are extremely useful in distinguishing 
different mechanisms of the transition. Further experimental studies are most desirable. 

There are a number of remaining problems. In the present analysis, we have not considered 
so far the second order processes with Ti — N intermediate states induced by the pion (meson) 
and/or quark exchanges. The short range part of weak transition potential for T,N NN was 
also computed in the same direct quark mechanism]^]. It is found that the mixing of this potential 
does not change the main feature of present DQ potential, though its contribution depend on the 
probability of the S mixture and is not negligible quantitatively. 

For hypernuclei other than the s-shell systems, we need a realistic calculation combined with 
the nuclear structure analysis. Further study is needed. 
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